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, Abstract 

>• 

' Numerical experiments on most gauge theories coupled with matter 

f**^ , failed to observe string-breaking effects while measuring Wilson loops only. 

f"^ ■ We show that, under rather mild assumptions, the overlap of the Wilson 

t-H ' loop operator with the broken-string state obeys a natural upper bound 

t-H , implying that the signal of string-breaking is in general too weak to be 

detected by the conventional updating algorithms. 

| In order to reduce the variance of the Wilson loops in 3-D Z2 gauge 

4-^ ■ Higgs model we use a new algorithm based on the Liischer-Weisz method 

I ' combined with a non-local cluster algorithm which allows to follow the 

I , decay of rectangular Wilson loops up to values of the order of 10 -24 . In 

this way a sharp signal of string breaking is found. 

^ : 

> ■ 1 Introduction 

V . 

. The confining force between a pair of static sources in pure gauge theories is 

mediated by a thin flux tube, or string, joining the two sources. When matter 



is added to this system the string becomes unstable at large separations and 
breaks when it reaches a certain length Rb to form pairs of matter particles. 
This breaking should produce the screening of the confining force between the 
sources and hence the flattening of the static potential. 

The lack of any sign of flattening in most systems, while measuring the static 
potential from Wilson loops only, came as a surprise E] • 

One suggestion to arise in literature is that the Wilson loop has a poor 
overlap with the true ground state, hence the basis of the operators has to 
be enlarged [3] in order to get a reliable estimate of the potential. Using this 
multichannel method it has been observed the breaking of the confining string in 
Higgs models g] in QCD 0E1 and in the SU(2) Yang Mills theory with adjoint 
sources [7]. Even if some cautionary observations have been raised about this 
method [5] , one is led to conclude that the difficulty in observing string breaking 
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with the Wilson loop seems to indicate nothing more than that it has a very 
small overlap with the broken-string state. 

This fact has been directly demonstrated in 2 + 1 d SU(2) YM theory with 
adjoint sources [H] where, using a variance reduction algorithm allowing to de- 
tect signals down to 1CP 40 , it has been clearly observed a rectangular Wilson 
loop W(R > Rb,T) changing sharply its slope as a function of T from that 
associated to the unbroken string (area-law decay) to that of the broken-string 
state (perimeter-law decay) at a distance much longer than the string breaking 
scale Here we will undertake a similar work for the 2 + 1 d I2 gauge- Higgs 
model and find a similar result. With the emergence of this long-distance effect 
a closely related question comes in: why the overlap of the Wilson loop with 
the broken-string state is so small?. 

Before trying to answer to this question, we should mention that do exist 
gauge systems where such overlap is much bigger. For instance, in the 2 + 1 d 
Z 2 gauge-Higgs ^H] an d in 3 + 1 d SU(2) gauge-Higgs models ^T] ^ nas been 
identified a region of the space of the coupling constants where the vacuum 
has a rather unusual property in that it has, besides the magnetic monopolc 
condensate characterising the confining "phase" , also a non- vanishing electric 
condensate, like in the dual Higgs phase. We argued ^U] that in this region the 
world sheet of the Wilson loop belongs to the so called tearing phase ^1 Ej , 
characterised by the formation of holes of arbitrary large size, reflecting pair 
creation. As a consequence, larger Wilson loops give rise to larger holes and the 
vacuum expectation value follows the perimeter-law, signalling string breaking. 
Actually, measuring Wilson loops in such a special region of the 2 + 1 d Z2 
gauge Higgs model, we observed ^H] a smooth transition from an area decay at 
relatively short distances to a perimeter decay at long distances. 

Hints of an appreciable overlap of large Wilson loops with the broken-string 
ground state have been also reported in 2+1 dimensional SU(2) gauge theory 
with two flavours |14j using highly anisotropic improved lattices and even in 
3+1 QCD 15,, measuring Wilson line correlators in Coulomb gauge with an 
improved action. 

In this paper we consider instead isotropic lattices with standard plaquette 
action in the confined "phase" . Our main goal is to understand why in these 
cases, for all coupled gauge systems studied thus far, it is so difficult to observe 
string breaking using Wilson loops only. 

It was argued ^1 that in these cases the world sheet associated to the 
confining string belongs to the normal phase, where the holes induced by dy- 
namical matter have a mean size which does not depend on that of the Wilson 
loop: larger loops give rise to the formation of more holes. As a consequence 
they decay with an area-law as in the quenched case. This suggests we assume 
the following Ansatz for the asymptotic behaviour of large, rectangular Wilson 
loops 1 

W(R, T) ~ c u exp[-2p(R + T) — aRT] + c b cxp[-2pt(i? + T)} (1) 

The first term describes the typical area-law decay of pure Yang-Mills theory 
with a string tension a . The second term is instead the contribution expected 
in the broken-string state; it decays with a perimeter law controlled by the mass 

x For sake of simplicity wc momentarily neglect the universal contribution of the quantum 
fluctuations of the flux tube. For an improved Ansatz see Ea. lf 61 below. 
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fjj of the so called static-light meson, the lowest bound state of a static source 
and a dynamical Higgs field or quark . 

In |13| it was even conjectured that the perimeter term could be zero (no 
overlap). Note however that any Wilson loop of finite size receives contributions 
not only from world sheet configurations typical of the normal phase, but also 
from those of the tearing phase [T§|, hence c& 7^ 0. 

It is easy to see [S] that for R < Rb and c b small enough the first term 
dominates over the whole range of T and no sign of string-breaking can be seen. 
When R > Rt, no matter how small Cb is, the above Ansatz implies that at long 
distances the broken string behaviour eventually prevails, since the first term 
drops off more rapidly than the second. We must of course also require that the 
latter should not dominate at distances shorter than the string breaking scale 
Rb, thus it is obvious that Cb cannot be too big. 

A closer look shows (Section |2J that the ratio Cb/c u is bounded from above 
by the following inequality 

log^ < aR f (2R b -R f ) , (2) 

where Rf is the minimal scale of string formation. Inserting this upper bound 
into Q one easily checks on the back of an envelope that as R ranges from, 
say, Rb to 2 Rb — Rf, the Wilson loop W(R, T) could deviate from the unbroken 
string behaviour only for T > 2 Rb — Rf , where its vacuum expectation value 
is in general exceedingly small (Section El . 

Thus such a bound gives us a simple explanation of the fact that it is so 
difficult to see a signal of string breaking while measuring Wilson loop only. 
It can also be used as a guide to search in the parameter space of the gauge 
theory a region where it is possible to evaluate the vacuum expectation value of 
rectangular Wilson loops with R > Rb and T > 2 Rb — Rf. 

We applied that procedure to the 2+1 dimensional Z 2 gauge-Higgs model 
(Section ISJ. Although such a system is perhaps the simplest example of a 
gauge theory coupled to matter in the fundamental representation, the current 
algorithms are not sufficiently accurate for our purpose. Thus we developed a 
new one, based on a version of the recent Liischer-Weisz variance reduction |17| 
for updating the gauge degrees of freedom, combined with a non-local cluster 
algorithm for the matter degrees of freedom ( Section 0}. 

The new algorithm allowed us to detect signals down to 10~ . We observed 
in this way rectangular Wilson loops W(R, T ~ 2 Rb) changing abruptly their 
slope as a function of R as described by the Ansatz Q}. It turns out that the 
breaking point nearly saturates the upper bound (Section EJ). 

2 The upper bound 

Consider a gauge system composed by a gauge field coupled to whatever kind of 
matter. The confining string between a pair of static sources should be unstable 
against breaking at large R, where dynamical matter particles can materialise 
to bind to the static sources, forming a pair of bound states called static-light 
mesons. 

Denoting by [i the mass of the lowest bound state, the static potential is 
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expected to approach the constant value 



Km V(R)=2fi , (3) 

R— too 

where 

V(R) = -]xm ^W(R,T) . (4) 

Although Eq refers to an asymptotic value, numerical simulations on these 
systems show that for all practical purposes we may safely assume 

V(R) ~ 2n , R > R b , (5) 

while in the range Rf < R < Rb, where the confining string is formed and is 
stable, the potential should have the typical confining form dictated by the first 
term of Eq. Q 

V(R) ~ 2p + aR , Rf < R < Rb ■ (6) 
Combining Eq.(|SJ) and JBJ) yields 

= ^ . (7) 
a 

Notice that the mass fj, and the perimeter term p are not UV finite because of 
the additive divergent self-energy contributions of the static sources. However, 
these divergences should cancel in their difference, hence Rb is a meaningful 
physical scale even in the continuum limit. 

In the range Rf < R < Rb, according to Eq.©, the first term of the Ansatz 
|T| dominates over the second in the limit T — > oo, hence it should also dominate 
for any finite T, because V(R) is less than 2 p in this range of R. Thus 

c u e - 2 K«+T)-<rAT > Cb e -2^R+T) 5 Rf < R < Rb , T > Rf . (8) 
With the help of Eq. Q , this inequality can eventually be recast into the form 
log— < a [R b (R + T) -RT] , R f < R < R b , T > R f . (9) 

Putting R = T = Rf yields the sought after bound For a graphical repre- 
sentation of the Ansatz and the subsequent bounds see Fig. ^ 

In the confining string picture, the inequality JSJ) tells us that a string of 
length Rf < R < Rb is stable and then it can propagate for any interval of T. 
A natural question comes to mind. What happens when we stretch the string 
beyond Rb? clearly it becomes unstable against breaking. Let T = T(i?) be 
the amount of "time" it takes for a string of length R to break, defined by the 
equality of the two terms of the Ansatz Q). Inserting there the upper bound 
@ we can find a lower bound for T as a function of R: 

Rb (R — 2 Rr) + Ri 
r{R) > R ^ f - , R> Rb ■ (10) 

T represents the minimal survival time of a string of length R before breaking. 
In order to maximise the signal one must minimise the area A(R) = RT(R) of 
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the Wilson loop. Looking at the solution R a of = one gets after a little 

algebra 

T b = T(R ) = 2R b — Rf = R . (11) 

It represents the most favourable distance for observing string breaking in the 
case in which the upper bound were saturated; but even in this limit case the 
signal is very weak, being proportional to e~ aTb . From a computational point of 
view it is very challenging to reach such length scales in the measure of W(T, R) 
even in the simplest models. This explains why it is so difficult to sec this effect. 

Thus far we have neglected the effects due to quantum string fluctuations. 
A more accurate Ansatz which accounts for these will be used while fitting the 
numerical data (see Ea. Qltjf) '). Had we used such a refined form in deriving the 
above bounds, we would have obtained much more involved formulae without 
modifying very much their numerical value. 

In deriving the previous inequalities it was assumed that the confining string 
is created by the Wilson loop. It should be clear that these formulae are not 
applicable to strings generated by different operators. Notice that different 
operators correspond to different boundary conditions of the string world sheet 
and these imply in turn different long distance behaviour. For instance, in the 
case of a coupled gauge theory at finite temperature, it is quite obvious how to 
modify the Ansatz (JJJ in order to describe the Polyakov line correlator (see e.g. 
Eq.(lO) of |13|). In particular the second term is just a constant at any fixed 
temperature; this implies that the overlap of the Polyakov line correlator with 
the broken string state is maximal above R b , in agreement with the fact that at 
finite temperature string breaking has been easily seen |20| . 



3 Z2 gauge-Higgs action and observables 

The action of a 2+1 dimensional Z2 gauge theory coupled to Ising matter in a 
cubic lattice A can be written as 

S(p a ,0r) = -fii ViUim -fc^Uu, (12) 

(ij) plaq. 

where both the link variable E/jy = Ue and the matter field ifi take values ±1 
and U a = l\ een U e . 

This model is self-dual: the Kramers- Wannier transformation maps the 
model into itself. Its partition function 

Z{Pg,Pi) = J [DU] \D<p] e -s(/Wj) (13) 

fulfils the functional equation 

Z(fia,0i) = (sinh2/3 G sinh2/3 / )i Ar Z(/3 7 , / 9 G ) (14) 

with /3 = — i log(tanh/3). 

The phase diagram of this model was studied long ago [T^] and revisited 
recently 1 19j. There is an unconfined region surrounded by lines of phase tran- 
sitions toward the Higgs phase and its dual confining phase. These lines are 
second order until they are near each other and the self-dual line, where first 
order transition occurs. Our simulations are of course in the confining phase. 
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We measured two kinds of observables. The Wilson loop associated to a 
closed path C of links £ is defined as usual 



W(C) = { J] U t ) = - f [DU] [Dvj] 1] Ut e~ s . (15) 

lec J lec 

We fitted the numerical data using a refinement of the Ansatz which takes 
into account the universal contribution of the bosonic string quantum fluctua- 
tions. Indeed previous numerical work on this subject strongly suggested |16j 
that the nature of the underlying asymptotic string should be the same as in 
the pure gauge model |21j . Thus we put 



W(R, T) ~ c u J ^ exp[-2p(R + T) - aRT] + c b cxp[-2^(R + T)] , (16) 
V w) 

where rj(r) is the Dedekind function 

OO rp 

V(r) = q 1/2i U(l q n ) , q = e 2 ^ , r = i- . (17) 

n— 1 

We also considered the gauge-invariant propagator G(T) of the static- light me- 
son. It can be constructed by coupling the product of link variables along a line 
of extent T to the matter fields ipo and tpT located at both ends 

G(T) = (<p Q ]lUt<PT) ■ (18) 

When this line is long enough the asymptotic form of its vacuum expectation 
value is 

G(T)~ce"^ T , (19) 

and is well suited to measure ix. We did not apply any kind of smearing to the 
measured operators. 



4 The algorithm 

The main idea underlying our algorithm relies on the consideration that the 
Liischer-Weisz procedure |17| . reducing the short wavelength fluctuations, is 
capable of an error reduction even if we cannot use their argument on the expo- 
nential decay of the temporal line correlators, because in presence of interacting 
matter the confining string breaks. 

We proceed as in the Liischer-Weisz algorithm and split the lattice in sub- 
lattices formed by temporal slices and evaluate Wilson loops via a stochastic 
estimate on these sub-lattices. To be specific, consider a rectangular R x T 
Wilson loop in the (x,t) plane. Let Uo(x,i) (Ui(x,t)) be the time-like (space- 
like) link variables. The two spatial sides of length R are associated to the 
operators L(0) = [] n= i R/ a Ui(na,Q) and L(T) = Yl n =i R/a Ui{na,T) and the 
operators associated to the two temporal sides are conveniently expressed in 
terms of pairs of time-like links with the same temporal coordinate t = na, 
which we call T(R,na) = Uo(0,na)Uo(R,na). We adhere to the notation of 
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Ref.|17|. but we need no matrix indices nor complex conjugation, of course. The 
vacuum expectation value defined in Eq. l|15|) can be rewritten as 

W(R, T) = (L(0) T(i?,?ia)L(T)) . (20) 

n=l,T/a 

Let us split the whole lattice in sub-lattices formed by temporal slices of arbi- 
trary thickness. As observed in ref.^7], owing to the locality of the action, each 
slice can be analysed independently of the surrounding medium, provided that 
the link variables of the boundary and the matter field configuration are held 
fixed. 

In the present case the expectation value of the operator T(i?) in a time-slice 
formed by to layers is defined by 

[T V {R, j;m}} = [T(R,ja)T(R, (j + l)a) . . . 

T(R, (j + m - l)o)] = ^—\ [DU] e -Wtf,d 

^sub J 

where the subscript sub refers to the variables belonging to the sub-lattice. 

In our numerical simulations these expectation values are estimated using 
a heat-bath method for the link variables, alternating in a suitable proportion 
with a non-local cluster algorithm for updating the matter field ip 

As in pure gauge case one gets identities like 

[T v {R, J]m }] [T v {R,(j + m);n}} = [T v {R,j;m + n}} , (22) 
which allow us to rewrite W(R,T) in the form 

W{R,T) = (L(0) [] [T v {RJ;m}]L(T)) . (23) 

j=l,T/m a 

where for simplicity it is assumed that T/a is a multiple of m. In our simulations 
we chose m = 2 and m = 3. 

An updating cycle is organised as follow: 

1. Generate field configurations of the whole lattice using the heath-bath and 
the non-local cluster updates as in |10| : 

2. Hold fixed both the spatial links along some suitably chosen spatial planes 
and the matter field configuration and update the gauge fields in the in- 
terior of the slices to estimate [T v {i?, j; to}]; 

3. Combine the stochastic estimates of the previous step with the side op- 
erators L of Eq. (|23|l and update the matter field configurations using the 
non-local cluster algorithm; 

4. Go to 1. 

The number n w of updates of the whole lattice (step 1), the number nt of updates 
inside each time-slice (step 2) and the number n v of updates of the matter field 
(step 3) are chosen in practice on the basis of a trial and error method in the 
optimisation process of the error reduction. For instance, in the to = 3 case one 
updating cycle was composed of nt = 10 2 gauge updates for each slice followed 
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Figure 1: Square Wilson loop versus R in a semilogarithmic plot. The parabolic 
line is the first term of the Ansatz |TSJ); the straight line represents the contri- 
bution of the broken-string state (the second term of the Ansatz) as determined 
by the fit to rectangular Wilson loops at T = 15 a and large R. 

by one non-local cluster update of the matter configuration (i.e. n v = 1) and 
n w = 10 5 gauge updates of the whole lattice. The total number cycles generated 
to extract our estimates was 4.0 x 10 4 . All the simulations for this work were 
performed on our cluster employing 20 processors for a total amount of 7000 
hours each. 

5 Results 

We first explored a wide region of the confining phase in order to see whether 
there exists a parameter range where the lower bound (|1C)|) is accessible to real- 
istic simulations. For this purpose we measured the vacuum expectation value 
of square Wilson loops W(T, T) and the static-light meson correlators G{T) in 
a wide range of T on a 40 3 cubic lattice in order to estimate the string tension 
a and the string breaking scale Rb trough Eq.(|7||. 

The square Wilson loop data were fitted to the first term of Eq. (jl6(l by pro- 
gressively eliminating the data of lower T until stable parameters were obtained. 
While doing this, the scale Rf of string formation was determined by picking 
out the value of T = Rf such that the exclusion (inclusion) of (W(Rf,Rf)) in 
the fit resulted in a good (bad) \ 2 test. In all the cases considered there was 
no ambiguity in the choice of Rf, being the corresponding variation of x 2 large 
enough. 

In this way we selected the point (3g = 0.650, (3i = 0.235 corresponding, in 
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Figure 2: Wilson loop data versus R, for T = 15 a. The steepest line is not a 
fit, but represents the first term of Ea. l|ro|) as determined by fitting the square 
Wilson loops. The slope of the other line is twice the mass of the static-light 
meson. The only free parameter is the intercept of this line. 

lattice spacing units a, to 

era 2 = 0.2054(21), fj,a = 0.8637(10), R f ~ 3a, R b ~ 7.5a, c„ = 1.19(5) . (24) 

Here we concentrated our computational efforts enlarging the lattice size to 50 3 
and measuring systematically all the rectangular Wilson loops with 4a < R < 
18a and T/a a multiple of 2 or 3. 

In Figlland FigEJwc show W(R, T) at T = 15 a and T = 14 a as a function 
of R. Notice the abrupt change in slope, signalling string breaking. The steepest 
line is not a best fit, but is drawn using the first term of Ea. i|16f) (unbroken- 
string term) where the parameters are those fitting the square Wilson loops. 
The slope of the other line is given by 2fi. The only parameter used to fit the 
data of Figgis Cb of Ea.Q16[l. We estimated c& = 43(5). This value is used in 
turn to account for the data of Fig|3J Inserting these values in Eq.® we see 
that the upper bound is nearly saturated and the string-breaking scale in Fig^ 
is correspondingly slightly bigger than T},. 

6 Conclusion 

The string breaking phenomenon in gauge theories coupled with matter is hardly 
detectable while measuring Wilson loops only. This reflects the poor overlap 
of the Wilson loop operator with the broken-string state. A simple way to 
represent such a behaviour is to assume that the vacuum expectation value of 
a rectangular Wilson loop W(R,T) is the sum of two terms (see Eq.Q), one 
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Figure 3: Wilson loop data versus R for T — 14 a. We did not fit any parameters 
but used those of the previous figure combined with the Ansatz JTSJ. 

decaying with an area law which prevails at intermediate scales and the other 
obeying a perimeter law which describes the expected asymptotic behaviour. It 
is worth noting that the logarithm of the former term is a hyperbola in the R, T 
plane, while the log of the latter is a straight line; they intersect at two points. 
One intersection represents the cross over to the broken-string state. The other 
intersection cannot be completely arbitrary: in order not to spoil the area law at 
intermediate distances one is forced to put an upper bound of the overlap to the 
broken-string state of the Wilson operator. Such an upper bound demonstrates a 
posteriori why it is so difficult to see string breaking in current simulations using 
Wilson loops only: even if the most favourable conditions were met, namely the 
saturation of our upper bound and the optimisation of the signal by measuring 
square Wilson loops, the minimal distance at which string-breaking is visible 
is Tf, = 2Rb — Rf (see Ea.(|ll|ll. where Rb is the string-breaking scale and Rf 
is the minimal scale of string formation. From a computational point of view 
Tf, is too large for the current updating algorithms. This explains why earlier 
studies on Wilson loops in gauge theories coupled to matter, which did not used 
a variance reduction method, failed to observe string-breaking. 

Adapting the Liischer-Weisz variance reduction method to the 3-D Z2 gauge 
Higgs model, which is perhaps the simplest gauge theory coupled to matter, we 
found a clean and beautiful signal of string breaking in the Wilson operator in 
a region where our upper bound is nearly saturated. 
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